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OPTIMAL RECURSIVE DIGITAL FILTERS FOR ACTIVE BENDING 

STABILIZATION 
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In the design of flight control systems for large flexible boosters, it is common 
practice to utilize active feedback control of the first lateral structural bending 
mode so as to suppress transients and reduce gust loading. Typically, active sta- 
bilization or phase stabilization is achieved by carefully shaping the loop transfer 
function in the frequency domain via the use of compensating filters combined 
with the frequency response characteristics of the nozzle/actuator system. In this 
paper we present a new approach for parametrizing and determining optimal low- 
order recursive linear digital filters so as to satisfy phase shaping constraints for 
bending and sloshing dynamics while simultaneously maximizing attenuation in 
other frequency bands of interest, e.g. near higher frequency parasitic structural 
modes. By parametrizing the filter directly in the z-plane with certain restrictions, 
the search space of candidate filter designs that satisfy the constraints is restricted 
to stable, minimum phase recursive low-pass filters with well-conditioned coeffi- 
cients. Combined with optimal output feedback blending from multiple rate gyros, 
the present approach enables rapid and robust parametrization of autopilot bending 
filters to attain flight control performance objectives. Numerical results are pre- 
sented that illustrate the application of the present technique to the development 
of rate gyro filters for an exploration-class multiengined space launch vehicle. 


1 INTRODUCTION 

The design of bending filters that robustly attenuate or active stabilize parasitic structural modes 
is of paramount importance in the development of a launch vehicle flight control system . 1,2 ’ 3 ' 4 The 
separation of the control input effector, typically a thrust vector control (TVC) system, and the trans- 
ducer element, a rate gyro, produces non-minimum phase parasitic elastic response with respect to 
thrust vector inputs. This structural response, in the case of large boosters, is highly coupled with 
other parasitic modes and appeal's at relatively low frequencies with respect to the control system 
bandwidth. Simple low-pass or notch filters commonly employed in aircraft feedback control appli- 
cations have insufficient performance to ensure robust stabilization of the structural modes, usually 
due to the penalty in open-loop phase that precludes stable control of the rigid body dynamics. Fur- 
thermore, the propellant sloshing modes may be in close proximity to the bending spectrum, further 
complicating the filter requirements. 

While rigorous controller design techniques such as H2 and H ^ 5 yield high-performance com- 
pensators, they do so at the expense of high controller order and limited design flexibility for a high- 
order plant, and do not guarantee classical stability margins . 6 In addition, since large boost vehicles 
exhibit uncertain nonminimum-phase elastic behavior over a broad spectrum of frequencies, the 
requisite model reduction of the plant or the resultant compensator can yield poor performance . 5 ' 7 
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Typical integrated vehicle linearized models used in control design arc of order 200 or more, and 
simultaneously consider all adverse control-structure interaction effects. 8 Although Hoc techniques 
have seen limited applications in launch vehicle flight control, 9 when algorithm complexity is con- 
sidered, experience has shown that classically tuned and highly optimized linear bending filters 
provide the best performance in this application. 1 - 10 

Methods of structural stabilization in launch vehicle control system design include active (phase) 
or passive (gain) stabilization. 3 In the case of active stabilization, a specific phase lag target is 
achieved by the bending filter in a region near the bending mode frequency. The open-loop peak 
response of the bending mode thus exceeds unity but does not encircle the —1 + j 0 point in the 
Nyquist plane. In the case of passive stabilization, the bending filter provides sufficient attenuation 
near the bending mode frequencies such that the open-loop gain of each passively stabilized mode 
does not exceed unity or some minimum value under all parameter perturbations. 

Active stabilization is easier to achieve for modes with high open-loop gain that are near the rigid 
body control frequency. The phase lag requirement of the filter is combined with the natural parasitic 
lag of the actuator dynamics, and is substantially less than that required to provide a high degree of 
stopband attenuation at the same frequency. Conversely, passive stabilization is often employed for 
high-frequency modes, where the uncertainty levels arc substantially higher, especially with respect 
to open-loop phase characteristics near the inertial coupling resonances or tail-wag-dog frequencies 
of the vectoring engines. Passive stabilization is easier to achieve for high-frequency modes due to 
the reduced phase penalty incurred by enforcing a stopband attenuation requirement at a frequency 
well above the rigid-body control frequency. 

In both cases, a numerical optimization technique is utilized to produce highly optimized in- 
finite impulse response (HR) filter designs that provide maximum robustness to structural mode 
uncertainty while minimizing phase lag at critical control system and sloshing mode frequencies. 
Various optimization and design methods have been employed for a variety of HR filter applica- 
tions, including spacecraft and launch vehicles. 1112,1314 In the following sections, the structure of 
a straightforward numerical optimization technique as well as its objectives and constraints will be 
detailed. This technique relies on direct pole-zero optimization in the complex ( z ) plane and is able 
to simultaneously satisfy phase constraints and attenuation objectives by employing the well-known 
sequential quadratic programming approach. 15 


2 FILTER OPTIMIZATION 


It is known that direct numerical optimization of transfer function coefficients is a very ill- 
conditioned optimization problem due to singularities in the gradient of the objective function. 
Given that the target filter is in discrete form, details of its structure can be presumed a priori 
so as to reduce the search space to a set of already-known feasible solutions. For example, we de- 
sire a stable, nonminimum phase SISO discrete transfer function. As such, all poles and zeros of 
the discretized system must lie strictly within the unit circle. Similarly to the technique presented 
in (Reference 16), we choose to parametrize each filter in zero-pole-gain form as 


_ fco U(z 2 - fijZ + ati) 
Y\{z 2 - fijZ + a.j) 


( 1 ) 


with i,j = 1 . . . k where k = n/2 and n is the filter order. Furthermore, [J'f — 4 cq < 0 and 
B'j — 4 ay < 0, so the poles and zeros of the discrete filter always appeal - in complex pairs. In this 
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form, we may introduce a transformation to polar coordinates in the z- plane. Let A* = r zl e ±J ° zt 
(A j = r p: j e ±jl>PJ ) be the roots of the i lh (j th ) zero (pole) polynomial in Eq. (1). We may then 
construct a parameter vector 

% — [ r zl 9 Z i Tp\ 0p\ . . . r z k 6 z i~ fpk 9pk ] (2) 

and impose hard limits on the parameters x that restrict the admissible solutions to yield reasonable 
filter designs having a low-pass characteristic. The gain parameter kg is always explicitly calculated 
to normalize the DC gain of the filter, H(z ) z= i = 1, where 

k 0 = H(z)/\H(z)\ z=1 . (3) 

In the most general sense, limit constraints on Xi impose sector constraints in the z- plane. For 
example, in order to yield stable, nonminimum phase filters, we require r zl , r p j < 1; in order to 
enhance robustness of the filter to truncation errors, the location of complex poles and zeros can be 
confined to a sector that is separated from the unit circle by some nonzero e. 16 - 17 An example of 
sector constraints are shown below in Figure 1 along with a typical filter pole-zero pattern. In this 
case, the zeros are confined to a sector closer to the boundary of the unit circle so as to ensure that 
all zeros are effective in shaping the response of the filter. 



Figure 1. Pole-zero sector boundary constraints 


The use of a transfer function for optimization of a general SISO linear system has the advantage 
of being a minimal realization. 18 If we consider the response of an n th order discrete SISO linear 
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system to be a nonlinear process as a function of its parameters x, 

H = H n (x,z), (4) 

and if the parameters x arc the transfer function coefficients, the 2 n parameters in x arc a minimal re- 
alization of //,, and completely span the objective space. Clearly, other discrete linear time-invariant 
realizations, such as the discrete-time state equations, arc not minimal realizations, but linear trans- 
formations can render the realization sparse and eliminate the redundant parameters. Examples of 
such transformations include conversion to phase variable canonical form or diagonalization. 

In the present application, it is neither necessary nor desired to span the entire objective space. 
By assuming that the class of discrete-time filters consisting of only complex pole and zero pairs 
will adequately satisfy the optimization objectives, one half of the search space can be eliminated 
via symmetry in the complex plane. It is also assumed that the number of zeros and the number of 
poles arc equal (the resultant transfer function is biproper). The addition of the hard limits on the 
parameters confines the search space to specific sectors in the complex plane, and the result is that 
only stable, minimum-phase filters with a prescribed characteristic are considered in the space of 
admissible solutions. The dimensionality of the search space remains the same (2 n) but singularities 
in the parametrization arc eliminated. 

An additional nonlinear constraint is valuable for optimization algorithm convergence and filter 
implementation stability. When performing finite differencing to compute gradients of the objec- 
tives, superposition of filter poles and zeros may yield undesirable “masking” of pole-zero pairs. In 
addition, it is known that sensitivity to coefficient truncation is approximately proportional to the 
proximity of the poles to the unit circle and to the system zeros. 17 In order to avoid these phenom- 
ena, a minimum norm distance constraint is enforced in the complex plane. Given n zeros \ and 
n poles Xj, taken pairwise, there are p = n\/ (2 (n — 2)!) unique pole-zero pairs. For relatively 
low-order filters as arc typical for bending filters, this is a modest number (p = {6, 28, 66} for 
n = {4, 8, 12}, respectively). Let P be the set of ah unique pole-zero indices. A supplementary 
constraint of the form 

c(x) = min ||Aj(x) — A_/(x)|| — 6 (5) 

i,jeP 

ensures that any adjacent pole-zero pair has the minimum separation distance 6, where 5 is a scalar 
parameter, and the constraint c(x) < 0. Due to symmetry, the constraint is applied only to the roots 
having positive complex parts. Since the index set P can be computed offline, determining c(x) is 
not a substantial computational burden even for orders above n = 12. 

With the present parametrization, the filter can be optimized as a general nonlinear goal attain- 
ment problem of the form 


min 7 subject to < 


F(x) — 7 w 
c(x) 


y %min 


< o 

< 0 


< x < x 


max 


( 6 ) 


where F{x) is a vector objective function, w is a goal weight vector, 7 is a scalar slack variable, and 
Xmin, 'J'rnax are the parameter limits. The nonlinear objective F(x) encodes filter design criteria, 
such as stopband attenuation, maximum passband ripple, phasing targets at particular frequencies, 
and so on, and the aforementioned constraint c(x) provides supplementary limits on the pole-zero 
separation. If all objectives cannot be simultaneously satisfied for a given filter order, it is often 
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pertinent to let those w, associated with the phase targets be zero (converting them to constraints), 
and weight the attenuation constraints according to their relative importance. Note that the number 
of objectives can be less than, equal to, or greater than the number of parameters, so the Jacobian 
V X F is not square in general. This type of goal attainment problem can be readily solved using 
finite differencing and general nonlinear multiobjective optimization methods, such as sequential 
quadratic programming (SQP). 1915 

3 DESIGN CRITERIA 

Based on conservative heritage design guidelines for boost vehicles, 3 filters are designed to at- 
tenuate structural flexibility, including all uncertainties, by a minimum of 6 dB below the critical 
gain in the case that bending modes arc passively stabilized. In the case of active stabilization, 
bending filters are designed to provide a minimum of 30 degrees of phase margin for phase stabi- 
lized bending modes when all uncertainties arc accounted for. The assumed viscoelastic damping 
in the structural model strongly influences the modal attenuation requirements imposed on the filter. 
All bending modes arc assumed to have a viscous damping ratio not exceeding one half percent 
since launch vehicle structural analysis and test has not demonstrated conclusive results supporting 
damping ratios above this level with flight-like boundary conditions. 

Many launch vehicle autopilots are of the simple proportional-derivative type, where the attitude 
errors arc a small-angle linearization of an error quaternion and the rate errors arc sensed based on a 
platform-mounted rate gyro. The open-loop frequency response of the controlled vehicle dynamics 
can therefore be expressed in the form 

— KpH att (juj )G a n (j/cj) T KdH rat (juj)G ra t {jui) (7) 

where I\ p , K,/ arc the proportional and derivative gains, respectively, and G a u and G ra t arc the 
open-loop system frequency response with respect to a thrust vector command input in terms of the 
attitude and rate errors. The frequency responses of the attitude and attitude-rate filter arc given by 
Hatt and H rat . 

Suppose a critical bending mode appeal's in the open loop at the frequency ui = loi>. If it is 
desired to phase-stabilize this mode, the total open-loop phase ZGorJj^) at the critical frequency 
uib should be near the critical phase 0i, = —3ir/2. Let the open-loop magnitude response at the 
critical frequency with H rat (juj ) = 1 be given by r = \G' OL (jub)\- To achieve a phase perturbation 
of the open-loop response without a substantial change in the open-loop magnitude at that frequency, 
the equality 

re j4>b = G 0 L{jub) (8) 

must be satisfied. For the purposes of determining the rate filter phase lag requirement, it will be 
assumed that the attitude filter is fixed. Equation (7) can be solved as 

H r at(jub) = (. K d G rat (juJb )) _1 (re j4>b - K p H att (ju b )G att {juJb )) (9) 

and the phase requirement for the filter at this frequency is simply the argument of the result. 

4 NUMERICAL EXAMPLE 

An example of the results of the present optimization method are shown below. The 8 th order 
digital filter design is implemented with a sampling period of 20 ms (50 Hz) and provides robust 
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Specification 

Rationale 

Frequency 

Value 

Phase 

Rigid-body phase margin 

0.20 Hz 

i t i> > —15 deg 

Phase 

First bending mode phase stability 

0.80 Hz 

-95 < <j> < -85 deg 

Phase 

First bending mode phase stability 

0.95 Hz 

-120 < <t> < -110 deg 

Phase 

Transition band phase shaping 

1.20 Hz 

(j> < —150 deg 

Attenuation 

Passband ripple 

0.10-0.60 Hz 

-1 < \H\ < ldB 

Attenuation 

Transition band peak 

0.60-2.00 Hz 

\H\ < ldB 

Attenuation 

High-frequency robust passive stabilization 

1.90-15.0 Hz 

\H\ < -26 dB 

Attenuation 

Notch high gain structural mode 

1.90-2.28 Hz 

\H\ < -32 dB 


Table 1. Example filter specifications 


phase stabilization of the first lateral bending mode of a large, multiengined exploration-class launch 
vehicle having a first lateral bending frequency near 0.9 Hz and a control system angular frequency 
of 0.2 Hz. The remaining structural attenuation is passive, and is achieved with strict attenuation 
requirements below the control system Nyquist frequency of 25 Hz. In addition, a particularly 
high-amplitude structural response near 2 Hz is suppressed by the addition of a supplementary 
attenuation constraint. The filter specifications and associated rationale arc shown in Table 1 . For 
implementation, the 8 th order optimal filter is decomposed into two 4 th order sections to avoid loss 
of precision on the target hardware. 

The response of the optimized filter is shown in Figure 2. The 8 th order filter is able to over- 
achieve the attenuation objectives while satisfying the phase shaping constraints that ensure rigid- 
body, slosh, and first bending mode stability. The resultant filter has additional robustness to uncer- 
tainty in the bending dynamics by increasing the stopband attenuation in the high-frequency portion 
of the spectrum. 
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Figure 2. Example filter response 


5 DISCUSSION 

The present parametrization technique offers a flexible design method which is compatible with 
standard software packages for nonlinear multiobjective optimization without overburdening the 
underlying numerical process. While the launch vehicle application drove specific requirements to 
shape phase and gain characteristics well outside of the capabilities of any existing analytic filter 
design methodology, this framework has seen successful applications in other aspects of aerospace 
vehicle signal processing. In addition to autopilot bending filter design, the software toolbox im- 
plementing this optimization method has been used to design downsampling filters to attenuate the 
resonant modes of sensor mechanical isolators, and to create linear approximations of noise pro- 
cesses based on their spectral descriptions. 
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